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Abstract 

In this paper, the Curie temperature of ferroelectric films is studied using spin- 
1/2 transverse Ising model with long-range interaction within the framework of the 
effective-field theory. The dependence of the Curie temperature on the thickness 
of the film, the surface interaction and the transverse field were investigated. It 
is assumed that the long-range interaction decays with the distance between the 
pseudo-spins as a power law. The dependence of the Curie temperature and the 
critical transverse field on the long-range exponent are obtained. 
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I. Introduction 

The pseudo-spin theory based on the Ising model with a transverse field (TIM) was 
first introduced by de Gennes [1] to describe the phase transition in order- disorder type 
ferroelectrics [2]. Under this model, it is assumed that a ferroelectric is composed of 
pseudo-spins with interactions, and this model had been applied to other systems suc- 
cessfully [3,4]. In most previous works [5-8], only the nearest neighbor interaction was 
considered. In ferromagnets where short-range interactions dominated, this simplification 
seems reasonable. But in ferroelectrics, where long-range interactions are thought to be 
dominant [9], the interactions between pseudo-spins over long distance should be taken 
into consideration. In previous works, the long-range coupling interaction has been con- 
sidered [10-12], it is found that the long-range interaction plays a vital role in ferroelectric 
structure. Theoretically, there are two methods to investigate properties of ferroelectrics, 
Ginzburg-Landau theory and the microscopic transverse Ising model. The long-range in- 
teraction in ferroelectric film has been studied by using TIM and the mean-field theory 
(MFT). However, the correlation between some spins is neglected in MFT. In order to 
study the critical behavior of ferroelectric film, some works have used TIM model within 
the frame work of the effective- field theory (EFT), however the long-range interaction 
was not considered [5-6]. It needs to use an approach, which is superior to the MFT. In 
this paper, we study the phase transition properties of ferroelectric film by use of EFT 
with correlation within the framework of TIM model, and the long-range interaction is 
considered in the model. 

The effective-field theory, which is based on the Ising spin identities and the differential 
operator technique, theoretically, is superior to the mean-field theory. The EFT has been 
used for the investigation of phase diagrams [5], and successfully applied to a variety of 
physical systems, such as magnetic thin films and superlattice [12,13], ferroelectric films 
[14] and ferroelectric superlattices [15]. In this paper, it is supposed that the interaction 
between the pseudo-spins decays as 1/r 6 [16], where 5 represents the decaying exponent. 
The case of 5 — > oo corresponds to the short-range interaction case, and the case of 5 — > 
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reduces to the infinite-range coupling one. According to our numerical results, we find 
that the long-range interaction has great influence on the phase transition temperature 
and the critical transverse field (i? c ). They both increase as the long-range interaction 
becomes stronger. 

II. Theory and Formulation 

We consider a film consists of N layers. Each layer is defined on the x-y plane and with 
pseudo-spin sites on a square lattice. As shown in Fig. 1, the 1-st and the N-th layers 
are the surface layers, in which the interaction constant is different from that in the inner 
layers. J s denotes the interaction constant of the pseudo-spins at the surface layer, and 
Jb that of the inner layer. We assume, for simplicity, that the interlayer interactions is Jb. 
In order to consider the long-range interaction in the ferroelectric thin film, we introduce 
the interaction constant J = Jij/rfj between pseudo-spins at i-th and j-th sites as in Ref. 
[16,17]. The system is described by the Ising Hamiltonian with a transverse field as 

ff — E^fsf-E^s* (i) 

where Sf and Sf are the components of the pseudo spin operator at site i, r^j is the 
distance between the i-th and j-th pseudo-spin, 5 is decaying exponent. 

According to the effective-field theory with correlation [18, 19]. The average value of 
pseudo spins in the i-th (i=l,2,3,...,N) layer can be given by 




II runs over all of the neighbors of site i, V = is the differential operator. The 
function F(x) is given by F(x) = ^tanh(Py), where (3 = -j^f, y = \J x 2 + i? 2 . 

The phase transition temperature is an important quantity in studying the critical 
properties of the ferroelectric film. Near the critical point, the polarization in each layer 
is very small, thus we can linearize Eq. (2). This leads to the matrix equation of the form 
Zr = 0. 
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For instance, we take N=10, and considering the 3-rd nearest- neighbor, then 



Rx = [cosh( J S V) + Ri sinh( J s V)] 4 [cosh( J h V) + R 2 sinh( J 6 V)] 

[cosh(^V) + R, sinh(^V)] 4 [cosh(^V) + R 2 sinh(^V)] 4 (3) 
a/2 a/2 a/2 a/2 

[cosh(-|y V) + it> 2 sinh( V)] 4 L =0 , 
i? 2 = [cosh( J s V) + i?i sinh( J S V)] [cosh( J 6 V) + i? 2 sinh( J fe V)] 4 

[cosh( J 6 V) + i? 3 sinh( J 6 V)] [cosh(^V) + i?i sinh(^ V)] 4 

a/2 a/2 

[cosh(-|,V) + it> 2 sinh( J,V)] 4 [cosh(J,V) + i? 3 sinh(-|_V)] 4 (4) 



[cosh( A_V) + Rl sinh (^V)] 4 [cosh(^V) + R 3 sinh(^V)] 4 F(x) 
Vo a/3 y3 a/3 



x=0 ' 



= [cosh(J s V) + R^ sinh( J s V)][cosh(J 6 V) + Ri sinh(J b V)] 4 

[cosh(J 6 V) + Ri +1 sinh(J b V)][cosh(^V) + sinht^V)] 4 

a/2 a/2 

[cosh(^V) + it> sinh(^V)] 4 [cosh(^V) + R i+1 sinh(^V)] 4 (5) 
v2 v2 v2 a/2 

[cosh(^V) + itU sinh(^V)] 4 [cosh( + it> +1 sinh(^ V)] 4 F(x)|_ , 

a/3 a/3 a/3 a/3 



i? 9 = [cosh(J s V) + R w sinh( J s V)][cosh(J b V) + R 9 sinh(J b V)] 

, [cosh( — ^tV) + Rg sinh , 
'a/2 ' V a/2° ' j 1 1 a/2 V a/2' 



cosh(J 6 V) + J R 8 sinh(J 6 V)][cosh(^V) + i? 10 sinh(^V)] 4 

a/2 a/2 

[cosh(^V) + i? 9 sinh(^V)] 4 [cosh(^V) + R 8 sinh(^V)] 4 (6) 
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-Rio — 



[cosh(^V) + R w sinh(^V)]>sh(^V) + R 8 sinh( A_V)] 4 F(x) \ x=0 , 
v3 v3 v3 v3 

[cosh( J S V) + i?i sinh( J s V)] 4 [cosh( J 6 V) + R 9 sinh( J 6 V)] 
[cosh(-^V) + R w sinh(-^V)] 4 [cosh(-^V) + it> 9 sinh(-^V)] 4 

[cosh(^ V) + it> 9 sinh(^ V)] 4 F(x) \ x=0 . 
v 3 v 3 



(7) 



Where i?j expresses the longitudinal polarization of a pseudo-spin in the i-th layer and 
the i takes the value from 1 to N. 

By linearizing the equations above, we obtain the matrix Z Ri = 0: 



where 
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) 



o, 



[4cosh(J fc V) cosh 4 ( A_ V ) cosh 4 (^ J V) cosh 4 (^V) sinh 3 (J s V) 



Jh 



V2' 



V2° 



V3° 



(8) 



Z\2 — 



+4cosh(J fe V) cosh 4 (J s V) cosh 4 ( A_ V ) cosh 4 ( A_ V ) sinh 3 (^V)]F(a;)| :!;= o, 

y/2 a/3 a/2 

[4cosh 4 (J s V) cosh 4 (^V) cosh 4 ( A_ V ) cosh 4 ( A_ V ) sinh(J 6 V) 



V2° 

+4cosh(J 6 V) cosh 4 (J s V) cosh 4 (^V) cosh 4 ( A_ V ) sinh 3 (^V) 



V2° 

A 
J, 



+4cosh(J fe V) cosh 4 (J s V) cosh 4 (— ^V) cosh 4 (^V) sinh 3 (^V)]F(a;)| :!;= o, 



V3° 

A 

V3° 

4 / Js 



3/ h 



V2' 
5 



V2° 



V3° 



x)\x=0- 

Lg the 



Z 21 = [cosh V) cosh 12 (^V) cosh 8 (^V)sinh(J 6 V) 

a/2 a/3 

+4cosh 6 (J 6 V) cosh 8 (^V) cosh 8 (^ V) sinh 3 (^ V) 

a/z a/3 a/2 

+4 cosh 6 ( J fe V) cosh 12 (^V) cosh 4 (^V) sinh 3 (^V)]F(a;)| :c= o, 

a/2 a/3 a/3 

Z 22 = [4cosh 2 (J 6 V)cosh 12 (^V)cosh 8 (^V)sinh 3 (J 6 V) 

a/2 a/3 

+4cosh 6 (J 6 V) cosh 8 (^V) cosh 8 (^V) sinh 3 (^V)]Fi 

a/2 a/3 a/2 

We can obtain the phase transition temperature of the ferroelectric film by takinj 

highest temperature from the solutions of the equation det(Z Ri ) = 0. 

III. Calculation and Discussion 

First, we putting Q j Jb = 0, 2.0, J s j Jb = I into the formulation in the previous section. 
Fig. 2 shows the T c versus n plot for different 5:5 = 1000 (in our calculation, it is 
equivalent to 5 — > oo approximately), 5 = 6.0, and 5 = 3.2. In Fig. 2, the solid line 
and dashed line correspond to Q j Jb = and Q j Jb = 2.0 respectively. For large 5(5 = 
1000), the long-range interaction is very weak, and it is corresponds to the short-range 
interaction. The smaller the value of 5, the stronger the long-range coupling interaction 
effect. From Fig. 2, we can see the dashed line of 5 = 1000 recovers the curve of 
transition temperature of a thin ferroelectric film versus the thickness n in the Zernike 
(Effective-field) approximation [5]. With the increase of the thickness, the phase transition 
temperature of the film approaches to the bulk one. 

As shown in Fig. 2, with the decrease of 5, which corresponds that the long-range 
interaction becomes stronger, the phase transition temperature increases. This is reason- 
able since strong long-range interactions represent that there are more other pseudo-spins 
around one pseudo-spin, which have the interaction on it. As a result, it can hardly flip 
freely. For the multilayer ferroelectric system, strong long-range interactions make the 
system harder to be disordered. Consequently, the transition temperature increases. 

Comparing Curve 2' in Fig. 2 in our paper with Curve 1 in Fig. 3(b) in Ref [11], we 
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find that they are very similar in appearance for the two curves. When film thickness 
reaches n = 10, the curve has almost tended to be saturate, i.e. T c tends to the bulk T c . 
It is clear that the MFT and the EFT obtain the same results. But the difference between 
Curve 1' in Fig. 2 in our paper and Curve 3 in Fig. 3(b) in Ref. [11] is obvious. The 
two curves both represent 8 = 3.2, but they are obtained by MFT and EFT, respectively. 
There is no T c = thickness in Curve 1' in Fig. 2 in our paper, while in Curve 3 in 
Fig. 3(b) in Ref [11], for a range of small film thickness in which the system has been 
disordered at zero temperature. 

It was very surprising that the Curie temperature and the polarization decrease with the 
increase of the long-range coupling interaction in Ref. [11]. Our result is that T c increases 
with the increase of the long-range interaction coupling as Fig. 2 shown. In physics, 
the interactions between the pseudo spins become stronger and stronger; therefore it is 
harder to be disordered for the system. Thus, only when the temperature is high enough, 
could the system become disordered. Therefore we believed that our results are more 
reasonable. 

Comparing the solid line with the dashed line in Fig. 2, we find the T c changes with 
different transverse field. Given the thickness n and the value of the decaying exponent 
5, the T c with a transverse field Q / Jb = 2.0 is lower than that without transverse field 
(See the dashed line in Fig. 2). i.e. the transverse field makes the T c of the system lower. 

The relationship between T c and the transverse field is shown in Fig. 3. We plot the 
T c versus i? with different 5. There exist a series of critical transverse field Q cs beyond 
which the system has already been disordered at K[5]. Q c decreases as the strength of 
the long-range interaction decreasing. The stronger the long-range interaction, the larger 
the critical transverse field Q c . 

Fig. 4 shows the Curie temperature as a function of surface coupling interaction J s . 
We take the thickness n=10. In Fig. 4a, we take different 5, which indicates the different 
range of the interactions. With the J s increasing, T c increases. At first, T c increases 
slowly. When J s reaches a certain point, T c increases linearly as the J s . We also find the 
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same result as in Fig. 2, i.e. smaller 5 results in larger T c . While in Fig. 4(b), taking 
5 = 6.0, we plot the T c versus J s in different transverse field f2/J b — 0.0, i?/Jf, = 2.0 and 
Q / Jb = 4.5. It is interesting that the curve of Q/Jb — 6.5 overlap the axis of J s when J s 
is smaller than a certain value. The reason is that: the transverse field Q suppress the 
orderness of the system. When transverse field is large enough, to a certain point, i.e. i? c , 
the phase transition point tends to zero for small J s , i.e. there is no ordered phase. While 
J s is larger than the value of J sc , the system can be ordered then. In ferroelectric film, it is 
of importance that we should take the long-range coupling interaction into consideration. 
The long-range interaction makes the transition temperature higher. 
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Captions of Figures 
Fig. 1. The schematic of a multilayer ferroelectric film 

Fig. 2. The temperature as a function of the film thickness n in the spin-1/2 Ising 
model with transverse fields ft / Jb = 0.0 (solid line) and ft/ J b = 2.0 (dashed line) 

Fig. 3. The phase diagram (T c versus ft plot) of the spin-1/2 10-layer system with 
a transverse field ft 

Fig. 4. The dependence of T c on the surface interaction constant J s . (a) for three 
selected 5 = 3.2, 6.0, 1000 with zero transverse field, (b) for three selected transverse 
fields: ft I J b = 0.0, ft / J b = 3.0, ft / J b = 6.5 with decaying exponent 5 = 6.0 
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